Capillary condensation in one-dimensional irregular confinement 
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A lattice-gas model with heterogeneity is developed for the description of fluid condensation in 
finite sized one-dimensional pores of arbitrary shape. An exact solution of the model is presented 
that reproduces the experimentally observed dependence of the amount of fluid adsorbed in the pore 
on external pressure. The proposed framework is viewed as a fundamental building block of the 
theory of capillary condensation necessary for reliable structural analysis of complex porous media 
from adsorption-desorption data. 
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Physical systems which consist of networks of pores, 
such as Vycor Silica aerogels porous rocks 
soil Q and others, have a wide spectrum of applications, 
ranging from molecular filters and catalysts [5| to fuel 
storage 0]. Capillary condensation is an important and 
interesting physical phenomenon occurring in many such 
systems. In recent years, a lot of experimental and theo- 
retical work has been undertaken in order to understand 
the link between the structure of porous materials and 
capillary condensation 0, 0, [!]• However, recent works 
suggest that understanding of capillary condensation is 
still poor The situation is critical since routinely 
used classical theories fail to explain experimental ob- 
servations 0] of capillary condensation even in simple 
onc-dimensional (ID) isolated pores. In particular, clas- 
sical theories have failed to explain the way in which the 
amount of adsorbed fluid varies as pressure increases and 
decreases, i.e. the appearance and form of hysteresis in 



sorption curves |10l. For example, according to clas- 
sical theories, pores closed at one end should not have 
hysteresis, which is in contradiction to existing experi- 
mental data [lo| . An accurate theory for capillary con- 
densation in ID pores is ultimately necessary to ensure 
that, for instance, this phenomenon can be used as a re- 
liable technique to probe the structure of generic porous 
media consisting of a network of interlinked ID pores [l^l ■ 

A possible theoretical explanation of these phenom- 
ena is that the appearance of hysteresis is caused by 
heterogeneity of the pores, which includes variable pore 
diameter, chemical heterogeneity in the pore walls and 
roughness of internal surfaces 3" 16 1- Numerical anal- 
yses such as lattice based mean-field theory 1171 and 
multi-scale molecular dynamics studies [l8|, |l9[ have 
given numerical support to this explanation. Addition- 
ally, these numerical approaches have revealed the occur- 
rence of avalanches (sudden jumps) in the amount of ad- 
sorbed fluid during condensation and evaporation. Such 
avalanches bare similarity to avalanches in magnetisation 
found in the zero-temperature random-field Ising model 
(zt-RFIM) [23, im . Employing this similarity, we map 



zt-RFIM to the lattice-gas model and demonstrate that: 
(i) a heterogeneous lattice-gas model is a minimally suf- 
ficient model to reproduce experimental observations in 
finite-sized pores; (ii) this model can be solved exactly 
analytically by a novel technique, with the solution being 
fully supported by numerical results; (iii) such an ana- 
lytical solution leads to simple physical explanations and 
interpretations of experimental results and predicts the 
shape of sorption curves for unexplored regimes. More 
specifically, we demonstrate that (i) the effects of the 
closed and open ends of the pores, considered impor- 
tant in the classical theories, are significantly reduced for 
large disorder in strengths of the interactions of the fluid 
with the pore walls; (ii) a positive-skewed distribution in 
strengths of such interactions can lead to sorption curves 
of the same form as those found experimentally [lol . 11 1 ; 
(iii) in cylindrical pores consisting of two parts of differ- 
ent diameter, condensation and evaporation in one part 
can induce condensation and evaporation in the rest of 
the pore for low disorder, but for high disorder, the two 
parts of the pore behave independently. 

The analyzed model is based on a standard lattice-gas 
model of capillary condensation 0, |2lj . In this model, 
the 3D space is split into cells, and each cell can be filled 
cither by matrix (the solid substrate), fluid or vapour. 
If the cell i is occupied by matrix (expressed by setting 
parameter rji = 0), then it cannot become occupied by 
fluid or vapour. It is assumed that the variables r]i are 
quenched for the whole system and thus the matrix state 
cannot change during the condensation. If cell i is not 
occupied by the matrix (rji = 1), then it can be occupied 
by either fluid (r^ = 1) or vapour (t^ = 0). The variables 
Ti can vary with the change in chemical potential, ^. 
The Hamiltonian which describes the lattice gas model 
is given by (Hj . 

"H = X! ^ X! TiTjTJ^TJj 
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FIG. 1. Diagram showing 2D sections of (a) a linear pore of 
N = 5 cells unoccupied by matrix, labelled i — 1, ... ,5 (white 
squares), and surrounded by matrix cells (grey squares) in 
the 3D system, (b) a hnear pore with open ends and (c) a 
linear pore of = 5 cells with one open and one closed end 
consisting of two sections represented by light- and dark-grey 
cells. The interaction of the fluid in cells 1 < i < A^i = 
3 and in cells 4 < i < A'' with the matrix represented by 
light-grey and dark-grey squares is described by {w^^), and 
(it)™'), respectively. For clarity, 2 relevant matrix cells per 
unoccupied cell, on top of and underneath the white cells, are 
not shown in this 2D diagram. 



where the summations run over all the cells in the system 
in the first term and over all nearest-neighbour pairs (ij) 
in the other terms. The fluid-fluid interaction parameter, 
> 0, is assumed to be the same for all pairs of cells. 
The matrix-fluid interaction strength between the matrix 
at cell j and fluid at the neighbouring cell i is described 
by the parameter wfj^. The values of u;™^ are consid- 
ered to be independently distributed quenched random 
variables with the probability density function, pi{wfj^), 
which can be cell dependent. The random distribution 
of this parameter has been studied previously in the con- 
text of chemical heterogeneity of the pore walls [22| , but 
below it is assumed to characterise all types of hetero- 
geneity. For concreteness, we focus on two forms of disor- 



describes the spin-spin interaction. The random fields hi 
at cell i are given by hi = J2j/iii^ - Vj)'^™^ I'^ + Vj'U'^ /^) 
with the sum running over all nearest neighbours of i. 
We study the zt-RFIM with cell-dependent random-field 
distributions. These distributions depend on the pore 
geometry chosen. Several forms of ID pore geometries 
embedded in a simple cubic lattice are analysed: linear 
pores (i) with closed ends (see Fig. [Tfa)), (ii) with open 
ends bounded by vapour (see Fig. [TJb)), and (iii) with 
one open end and one closed end. In the last case, we 
consider two distinct sections of lengths A''i and N — Ni 
characterised by different values of mean matrix-fluid in- 
teraction (w™^) and (w™^), respectively (see Fig. [Ijc)), 
intended to model ink-bottle and funnel geometries, with 
weaker matrix-fluid interaction representing a larger di- 
ameter. The random field, hi, at cell i is characterised by 
a distribution, ph.[hi), which depends on the numbers, 
= Ej7j(1 - V3) and 71^ = J2j/tVj, of neighbour- 
ing cells which are occupied and unoccupied by matrix, 
respectively. When wfj^ is normally distributed, the ran- 
dom field at cell i is ciistributed according to pinihi) = 
AA((/i,)„A,,.2) with (h,), = «/2)(u;-f). + «/4)«;* 
and A/j.^ = (n"y2)A|. When wfj^ follows a correlated 
exponential distribution, the random field at cell i is dis- 
tributed according to Phi(hi) — 6(?/)A^.^ exp(— A^.^y) 
(where y = hi — {hi)i + A/j.), with the mean {hi)i being 
the same as for the normal distribution and standard de- 
viation Afii = n™Ai/2. According to the dynamical rules 
of the zt-RFIM, for adsorption (desorption) each cell is 
initially empty (occupied), Si = —1 {si = +1), and can 
become occupied by fluid (empty), = +1 (sj = — 1), 



H 



IS positive. 



once its local field, fi = hi 
fi > (negative, < 0). 

In order to find sorption curves for a ID pore, we de- 
velop a novel method based on the use of the gener- 
ating function formalism [23|, |2J]. Within this formal- 



ism, the volume of fluid in the finite pore is a random 
value and can be expressed by a generating function, 

N 



der in wf"^, representing heterogeneity on different scales ^ix) = J2n=o Pin)x'^, where P{n) is the probability that 



ranging from local fluctuations at a single point on the 
pore wall to a variable diameter of the pore: (i) a nor- 
mal distribution, pi{w^^) = A/'((w™^)i, A^) and (h) an 
exponential distribution with correlations ensuring that 
all wfj^ are the same for the same cell i, and distributed 
according to, pi{'wfj^) = &{x)A~^exp{—A~^x) (where 



"^)i) -I- Ai and Q{x) is the Heaviside step 



function), with mean {w'^^)i and variance Af. Thermal 
fluctuations are ignored, under the assumption that the 
temperature is much smaller than all typical energy scales 
in the problem, T < w^, {w"'^)i. 

The Hamiltonian described by Eq. ([1]) can be mapped 
onto the Hamiltonian for the zt-RFIM [2lj, i.e. Ti — 
'^'^J2{ij) ^i^j ^ J2i^i^i ~ HJ2iSi, where the variable 
Si — {2Ti — Vj'qi {si — ±1) represents a spin state, H = 
p/2 refers to an external magnetic field, and J = w^/A 



'l2iLi ''"i cells in the finite chain are occupied by fiuid. 
The mean volume of fluid is given by 

{V) ^ Vo{n) ^ d.,G{l) , (2) 

with volume of a cell set to Vb = 1, and the variance is, 

Var[y] = a,,G(l) + 9,G(1) - [a,G(l)]' , (3) 

where dxG{l) and dxxG{l) refer to the first and second 
derivatives of G{x) with respect to x evaluated at a; = 1. 

Based on the recursive relaxation scheme used in 
Ref. (25! , an exact expression for the generating function 
G{x) can be derived (see [26| for technical details). 



Gix) = [Aix)y MN-iix)MN-2ix)...M2ix) [B(x)] 



(4) 
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where, 

[A{x)f ^ (xp^ 1 + (1 -p^i), a;p^o, 1 -p'n,o) 
M,{x) = 

xpi.i xpifi 

x{Pt,2 ~ Pt,l) + {'^ - P1.2) x{pis-Pi,o) l-Pi,i 



(a) Normal, A = 0.1 

Ir 



) Exp, A = 0. 1 



(c) Both, A = 0.25 

Ir 



Here, the values of p,^ „ = p[^^ (i = 1, iV) and p^,™ = p^^l^ 
{1 < i < N) are the probabilities that cell i, with m 
occupied neighbouring cells, has a positive local field, 
which arc given by p[% = lZ-H-Ji2m-k) Ph,ih)ih. 
Eq. is the main analytical result of our analysis al- 
lowing exact evaluation of dxG{l) and dxxG{l). The 
mean and variance of the volume of fluid in the pore 
can be found for both adsorption and desorption regimes 
using Eqs. ([2]) and ([3]) along with the derivatives of 
G{x). Technically, the derivatives of G{x) can be cal- 
culated by numerical iteration, i.e. the derivatives of 
Mi{x)Mi-i{x) . . . M2{x)B{x) can be found in terms of 
the derivatives of Mi-i{x) . . . M2{x)B{x). 

The analysis presented below is based on the numer- 
ical solution for G{x) given by Eq. ^ and supported 
by Monte-Carlo simulations of condensation in a pore 
within the framework of the lattice-gas model. During 
each simulation, the value of fj, was fixed, and hystere- 
sis plots were obtained by running separate simulations 
for each value of The fluid occupation number of the 
matrix-free cells was changed following single spin-flip 
zero-temperature Glauber dynamics 



20, 27 



The mean 

and variance of V were calculated by averaging over 10^ 
realisations of the disorder in random-fields for a fixed 
matrix structure and values of the parameters. Figs. [2] 
and [3] show that the agreement between analytical cal- 
culations (lines) and numerical simulations (symbols) is 
excellent. 

Our first result is that, for low disorder in w™^, ad- 
sorption in finite open-ended pores occurs in a single 
avalanche. This is evident from the fact that the variance 
reaches almost the maximum value Var[T^]i„ax = 0.25 (see 
right dashed peak in the lower panel of Figs.[2l^a)) imply- 
ing that the system is either fully occupied or fully un- 
occupied and there are no intermediate stable states (no 
pinning points). In contrast, for pores with closed ends, 
adsorption occurs at a lower value of /x (compare right 
hand solid and dashed curves in upper panel of Figs.[2l^a)) 
with the propagation of a meniscus between several pin- 
ning points in a number of separate avalanches, indicated 
by a lower Var[y] (see right solid peaks in the lower panel 
of Figs. [2ja)). The reason for this is that one of the 
cells at the ends of the closed-ended pores is occupied 
first with high probability, thus nucleating the adsorp- 
tion process in the form of meniscus propagation. For 
the open-ended pore, no cell can become occupied until 



>0.5 




FIG. 2. The mean (V) (upper panel) and variance Var[V] 
(lower panel) of occupied volume of linear pores with length 
A'^ — 100 is plotted vs fi. The arrows show the direction in 
which jj, changes for adsorption and desorption curves. The 
solid (dashed and dot-dashed) curves correspond to pores 
with closed (open) ends. Different columns refer to differ- 
ent types and degrees of disorder in wfj^ with the same mean 



value (w""') 



1.0 and the same w 



ff 



1.0: (a) p(w;'"*) 



normal with width A = 0.1; (b) correlated exponential dis- 
tribution with dispersion A — 0.1 (c) normal (dashed) and 
correlated exponential (dot-dashed) with width A = 0.25. 
Symbols refer to the numerical results. 



a much higher value of fi is reached, so that there are no 
pinning points at which the propagation of meniscus can 
stop and the avalanche, initiated at any cell in the pore, 
proceeds along the whole pore in one step. For desorp- 
tion, the situation is reversed to that for adsorption, with 
nucleation (evaporation) occurring at the end-point cells 
i = 1, N hi the open-ended pores first. At high disorder, 
both adsorption and desorption in all geometries shown 
in Fig. [1] occur mainly via many small avalanches (indi- 
cated by a low variance in the lower panels of Fig. EJc)). 
The sorption curves for open and closed ended pores are 
identical, agreeing with numerical results found for 
pores with both ends open or one end open. This is 
because there are some cells with high which nucle- 
ate the adsorption and pin the desorption and some cells 
with low which pin the adsorption and nucleate the 
desorption. 

Our second result is that the correlated exponential 
distribution of disorder in w™^ can result in a so-called 
H2-type hysteresis [l^ (see upper panel of Fig. Wljo) and 
dot-dashed lines in upper panel of (c)). This effect can 
be attributed to asymmetry in the distribution of ran- 
dom fields hi , because a symmetric distribution will lead 
to a symmetric hysteresis loop, due to the symmetry of 
the Hamiltonian of the zt-RFIM. The asymmetry in the 



distribution of /i^ = ^ 



mf 



can be achieved by a cor- 



related asymmetric distribution of wfj . Indeed, if wfj 
are uncorrelated for a given value of i then the central 
limit theorem ensures that their sum is approximately 
distributed according to a normal (symmetric) distribu- 
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(a) Normal, A = 0.1 




(b)Exp, A = 0.1 
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(c) Exp, A = 0.25 
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FIG. 3. iy) (upper panel) and Var[V^] (lower panel) for the 
linear pores shown in Fig. [TJc) of length A'^ = 100 vs ^. The 
curves of different styles refer to different geometries, i.e. ink- 
bottle (solid lines), with A''i = 50, liii = 1.2, W2 ~ 1-0 and 
funnel (dashed lines), with A'^i = 50, wi — 1.0, W2 = 1.2. The 
light dotted lines in the upper panels of (a)-(c) correspond to 
the mean occupied volume of two separate open-ended pores 
(shape (b) in Fig. [T]) of length JV = 50 one of which having 
(lu™') — 1.0 for all cells and the other having (wfj^) = 1.2 for 
all cells. In (c), the solid, dashed and dotted lines coincide on 
the scale of the graph. In all cases, = 1.0. Each column 
represents a different degree or form of disorder: (a) normally 
distributed, A = 0.1 (b) correlated exponentially distributed, 
A = 0.1 and (c) correlated exponentially distributed, A — 
0.25. Arrows show the direction of change of jj, for adsorption 
and desorption and symbols refer to numerical data. 



tion. Therefore correlations in wf-^ play a signifieant role 
in aehieving a skewed distribution of hi and H2-type hys- 
teresis, the effect being maximal when they are fully cor- 
related, i.e. when all are equal for a given i. This 
implies that H2-type hysteresis in heterogeneous pores 
arises due to variations in pore diameter (represented by 
correlated disorder in wf^^) rather than due to individ- 
ual defects, in agreement with previous numerical stud- 
ies 



17[. On the other hand, symmetric (normal) disorder 
in local fields (i.e. uncorrelated disorder in wfj^) can only 
cause a parallel sided Hl-type hysteresis, and may repre- 
sent defects on small scales (Fig. [2l^a)). 

Our third result refers to an explanation of the ex- 
perimental observations of condensation in pores with 
ink-bottle and funnel geometries. In pores consisting of 
two connected cylindrical sections of different diameters 
and a rough surfaces, it has been observed [llj that the 
condensation of fluid in the two parts of the pore occurs 
independently and is unaffected by which end of the pore 
is open. We have found that a similar effect occurs for the 
ink-bottle shape (see Fig. [TJc)) where the sections with 
different values w™^ and w™^ represent sections of differ- 
ent diameter used in experiments. Indeed, in Fig. [3jc) 
the adsorption desorption curves are shown for large, 
correlated disorder in w™^. The resulting curves are 
identical to those found by adding together the adsorp- 
tion desorption curves for the two isolated halves of the 



ink-bottle pore. For completeness, we also present the 
adsorption-desorption curves for small degrees of uncor- 
related (Fig.[n{a)) and asymmetric correlated (Fig. [3Jb)) 



disorder in 



In this case, there are significant in- 



teractions between the two sections, meaning that ad- 
sorption and desorption occur in the wider section at a 
higher value of /x in the funnel geometry (dashed lines in 
Fig.[3lja)-(b)) than in the ink-bottle geometry (solid lines 
in Fig. Ha)- (b)). 

To conclude, a heterogeneous lattice-gas model has 
been proposed to describe fluid condensation in ID pores 
of different shapes and rough surfaces. Heterogeneity is 
the key and sufficient feature of the model which allows it 
to reproduce the main experimental findings. In physical 
systems exhibiting a more complex topology, e.g. a 3D 
maze-like network of ID channels, an analytical solution 
may be nontrivial, however numerical simulations within 
this model can be performed straightforwardly for porous 
media of arbitrary topology. 

TPH acknowledges the support of EPSRC. 
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In this supplementary note, we derive an exact expression for the mean and variance of the volume 
occupied by fluid in a ID pore consisting of N cells. This is done by establishing a recursion relation 
for the generating function, 

N 

G(x) = ^P(nK, (1) 

n=0 

for the probability, P{n), that n of the cells in the pore are occupied by fluid. The abelian property 
of the zt-RFIM means that the final state of the system is independent of the order in which the cells 
are occupied by fluid (the order spins flip) as long as the system ends up in a stable or metastable 
state Q]- Therefore, we use particular dynamical rules for occupancy of cells that are convenient from 
a mathematical viewpoint. Assume that the relaxation of the system into a metastable state takes place 
in a series of time-steps, t, 1 < t < N . Initially, during time-step t = 1, all the cells with i > 2 are 
artificially prevented from being occupied, while the first cell i = 1 is allowed to change its state. Cell 
i = 1 can either change state from unoccupied to occupied if the local field is positive fi > 0, or remain 
unoccupied if the local field is negative. This process with two possible outcomes is called relaxation of 
cell 1. Next, during time-step t = 2, we allow cell z = 2 to relax, while cells z > 3 are still held in the 
unoccupied state, and cell i = 2 can become occupied if /2 > 0. If cell i = 2 does become occupied, 
then the local field /i at cell i = 1 will increase. This can cause cell i = 1 to become occupied if it was 
not occupied already, i.e. an avalanche can pass from cell i = 2 to cell i = 1 during time-step t = 2. 
Similarly, we then allow the next cell in the chain, i.e. i = 3, to relax and if it becomes occupied an 
avalanche can pass back along the pore towards cell i = 1 if those cells with i < 3 were not occupied. 
This method is recursively applied until all the cells in the system are relaxed, after time-step t = N. 

Let us consider the cell i = N at the end of the pore, which is the last cell to be allowed to relax 
in the above procedure. At the start of time-step t = N, the neighbouring pore cell i = N — 1 can be 
occupied or unoccupied, i.e. sjy-i = ±1- If the neighbouring cell i = A — 1 is occupied {sn-i = +1), 
and the random field at cell i = N is > — J — H, then the local field /at > and cell i = N will 
become occupied. This occurs with cell-dependent probability p'j^ ^ where, 

/•CO 

Pi,m = / pRF,iihi)dhi , (2) 

Jhi=-H-J{2m-1) 

with m being the number of occupied neighbours of cell i = N, i.e. m = 1 in this case. If, however, 
the neighbouring cell z = A — 1 is unoccupied {sm-i = — l), then the local random field must be above 
a higher threshold for cell A to become occupied, h^' > J — H, which occurs with probability, p'jy^, 
given by Eq. ([2]) with m = 0. In this case, the field at cell i = N — 1 will increase, and an avalanche can 
propagate back along the pore. 

The probability that there are n occupied cells in the lattice at the end of the relaxation process can 
therefore be written as, 

P(n) = Pa{N - 1, n - + P^(A - 1, n)(l - p'^^^) 

+ Ps{N-l,n-l)p'^^o + Pc{N-l,n){l-p'^^o) , (3) 

in terms of the probabilities, PA{i,n'), PB{i,n') and Pc{i,n'), which can be recursively determined. 
The quantity, 

PA{i, n) =Prob[ni(t = i) = n n Si{t = i) = +1] , (4) 

is the probability that at the end of time-step t = i, there are ni{t = i) = n' occupied cells j in the 
range 1 < j < i (i.e. = J2]=ii^j + l)/^) and that cell i is occupied, Si{t = i) = +1. The value, 



PB{i,n') = Prob[si(t = i) = -lr\ni{t = A) = n'\si+i{t = A) = +1] , 



(5) 
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is the probability that ceh i is unoccupied at end of time-step t = i, Si{t = i) = —1, and at the end of 
the relaxation process (end of time-step t = N) there are rii = n' occupied cells in the range 1 < j < i, 
given that cell i + 1 becomes occupied, Si+i(t = A^) = +1, during some time-step t', i < t' < N (which 
causes avalanches to pass back along the chain towards cell 1, changing the occupation number rii). The 
third quantity, 

Pfj{i, n) =Prob[ni(t = i) = n' r\ Si{t = i) = -I] , (6) 

is the probability that at the end of time-step t = i there are nj(t = i) = n' occupied cells in the range 
1 < i < ^ and that cell i is unoccupied at this time-step, Si{t = i) = —1. Using Eq. (I3|), the generating 
function G{x) for the total number of occupied cells in the pore can be written in terms of the generating 
functions Ai{x), Bi{x) and Ci{x) for the corresponding probabilities defined by Eqs. dH-dS]) as, 

G{x) = [Vat,! + (1 -p'n,i)] An-i{x) + xp'^^qBn-i{x) + (1 - p'^^q)Cn-i{x) . (7) 

where the generating function Ai{x) is defined as Ai{x) = X]n=o "')^"' while Bi{x) and Ci{x) are 
defined according to the same relation with A replaced by B and C, respectively. 

Expressions for the generating functions Ai{x), Bi{x) and Ci{x) for i > 1 can be found recursively 
in the following way. If cell z — 1 is occupied at the start of time-step t = i and the local field at cell i 
is positive, i.e. fi = hi + H > 0, then cell i will become occupied during time-step t = i. This occurs 
with probability pi^i, where 

POO 

Pi,m. = / ph,{h)dh . (8) 

Jh=-H-J{2m-2) 

If cell i — 1 is unoccupied then the random field has to be at a higher threshold hi > —H + 2J in order for 
cell i to become occupied during time-step i. The random field will be above this higher threshold with 
probability pi^, given by Eq. ([8]). If cell i does become occupied during time-step t = i and Sj_i = — 1 
at this time, then an avalanche might pass back along the chain towards cell i = 1. This gives the 
expressions, 

PA{i, n') = PA{i -l,n' - + Peii -l,n' - l>i,o , 
Pcii,n) = PAii-l,n)il-pi,i) + Pcii-l,n')il-p,^o) ■ 

(9) 

If cell i does not become occupied at time step t = i, but cell i + 1 becomes occupied at some later 
time-step t' {i < t' < N) then cell i can also become occupied during time-step t' . The probability 
of cell i becoming occupied in this way at time step t' depends upon whether cell i — 1 is occupied or 
not. In fact, if cell z — 1 is occupied then cell i becomes occupied when cell i + 1 becomes occupied 
only if the random field at cell i is in the range —H < hi < —H + 2J, which occurs with probability 
Pi,2 — Pi,i- On the other hand, if cell i — 1 is unoccupied then the random field must be in the range 
—H — 2J < hi < —H in order for cell i to become occupied when cell i + 1 becomes occupied, which 
occurs with probability pi^i —pifi- In the case that cell z — 1 is unoccupied, an avalanche of spin flips can 
propagate back down the chain from cell i towards cell 1 during time step t'. This gives the expression, 

Psii, n') = PA{i - 1, n')(l - ^,2) + PA{i - 1, n' - l)(p,,2 - Pi,i) 

+ Pc{i - l,n')(l + Pij(i - l,n' - l)(p,,i - Pi,o) • (10) 

Eqs. and (jlOp lead to the following recursive relations for the generating functions, Ai{x), Bi[x) and 
Q(x), 

Ai{x) = X [Ai^i{x)pi^i + B,^i{x)pifl] 

Bi{x) = x{pi^2 - Pi,i)Ai^i{x) + (1 - pi^2)Ai-i{x) + x{pi^i - pifi)Bi^i{x) + (1 - pi,i)Ci_i(x) 
Ci{x) = (1 - pi,i)Ai^i{x) + (1 -pi,o)C,„i(x) , (11) 
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valid for i > 1. 

The boundary values of P^(l,n') and Pc'(l,n') can be found using the following relations, 

PA{l,n') = 5n',lPi,o > 

= V,o(l-pi,o) , (12) 

where p'l q is the probability that cell 1 has a positive field (and thus becomes occupied) at the first time 
step (when all other cells are unoccupied). The value of Psi^^n') is given by the relation, 

PB(l,n') = (1 -p;,i)5„,,o + (P'i,i -p'i,o)V,i , (13) 

where p'l ^ — p'iq is the probability that cell 1 has a negative local field during time-step t = 1, but the 
field becomes positive when cell 2 becomes occupied, and 1 — p'n is the probability that cell 1 still has 
a negative local field after cell 2 becomes occupied. Eqs. (fT2]l and (fT3|l result in the following expression 
for the boundary generating functions, 

Ai{x) = xp{q , 

Bi{x) = a;(p; 1 - p;_o) + 1 > 

Ci(x) = i-p;,o- (14) 

The generating function given by Eqs. ([7]), ()lip and ()14p can be written as a matrix equation, 

G{x) = [A(x)]^ Mn-i{x)Mn-2{x) . . . M2{x) [B{x)] , (15) 



where. 



A{x)]^ = (Vat,! + (1 -p^_i), xp^o, 1 



' 



2;pi,i xpifl 

Mi(x) = I x{pi^2 - Pi,i) + (1 - Pi,2) x{pi^i - pifl) 1 - pi^i I , for 1 < i < iV , 
1 - Pi,i 1 - pi,oy 

[S(x)]T = (Vi,o, ^(Pi,i - Pi,o) + 1 - 1 - ^'i.o) • (16) 
Eqs. ([15]) and (fT6|) are the expressions quoted in the main text. 
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